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We investigate the tunneling properties of collective excitations through a domain wall in the 
ferromagnetic phase of a spin-1 spinor Bose-Einstein condensate. Within the mean-field theory at 
T = 0, we show that the transverse spin wave undergoes perfect reflection in the low-energy limit. 
This reflection property differs considerably from that of a domain wall in a Heisenberg ferromagnet 
where spin-wave excitations exhibit perfect transmission at arbitrary energy. When the Bogoliubov 
mode is scattered from this domain wall soliton, the transmission and reflection coefficients exhibit 
pronounced non-monotonicity. In particular, we find perfect reflection of the Bogoliubov mode 
at energies where bound states appear. This is in stark contrast to the perfect transmission of the 
Bogoliubov mode with arbitrary energy through a dark soliton in a scalar Bose-Einstein condensate. 
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I. INTRODUCTION 

Bloch and Necl walls in ferromagnetic materials are do- 
main walls with well-known structures in condensed mat- 
ter, where spin configurations are twisted in a magnetic 
domain region with an angular displacement of 180° [J. 
It is desirable to use domain wall motions in ferromagnets 
in the aim of their application to memory devices |2h13| . 
For example, current-driven domain wall motion associ- 
ated with spin/momentum transfer from electrons to the 
domain wall has been theoretically discussed p} and has 
been experimentally observed [|| . 

Domain walls arc also found in two-component Bose- 
Einstein condensates and spin-1 Bose-Einstein conden- 
sates. For spin-1 Bose-Einstein condensates, spin do- 
main structures [3 [l5| , tunneling from metastable spin 
domains 1161 1 . and the formation of spin textures and 
domains [171 Il8j have been experimentally investigated. 
An earlier theoretical paper [l9[ reported that ferromag- 
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FIG. 1: (color online) One-dimensional tunneling of excita- 
tions through a domain wall of a Bose-Einstein condensate. 
For a ferromagnetic spin-1 Bose-Einstein condensate, A and 
B represent different hyperfine states of the condensate (A 
corresponds to the S z = — 1 state and B to 5 2 = +1). A 
right-moving incident excitation is scattered from the domain 
wall. 
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netic domain walls in a spin-1 Bose-Einstein condensate 
spread out over the entire spatial region of the system. 
However, a domain wall may exhibit solitary behavior 
with a size comparable with the healing length of the 
Bose-Einstein condensate, which is much smaller than 
the system size. In addition, a uniform particle density 
around the domain wall was reported in Ref . fl9j , which 
is analogous to that of the Heisenberg model, where the 
magnitude of the spin vector is spatially uniform. How- 
ever, gaseous ferromagnetic spin-1 Bose-Einstein conden- 
sates may be free from such a constraint, so that their 
particle/spin density may become non- uniform in the do- 
main wall region. A search for such a domain is antic- 
ipated to reveal a solitary domain wall in spin-1 Bose- 
Einstein condensates. 

Furthermore, transport phenomena in gaseous ferro- 
magnetic Bose-Einstein condensates are much less under- 
stood than those in solid materials. Scattering between 
low-lying modes and a domain wall in a Bose-Einstein 
condensate will be interesting, if it exhibits different char- 
acters from conventional ferromagnets. In contrast to the 
Heisenberg model case, the spin density vector amplitude 
may be spatially non-uniform in spinor Bose-Einstein 
condensates. This character is expected to provide in- 
teresting insights into ferromagnets and spinor Bose- 
Einstein condensates. In addition, since these scattering 
problems can be expressed in terms of Nambu-Goldstone 
modes and topological excitations, this concept is not 
specific to cold atomic gases, but can also be applicable 
to other symmetry-broken states. 

We investigate the tunneling properties of collective 
excitations through a ferromagnetic domain wall in a 
spin-1 Bose-Einstein condensate. We start by reconsid- 
ering the domain structure of the ferromagnetic Bose- 
Einstein condensate with a given boundary condition. 
For simplicity, we consider a three-dimensional system 
with a planar domain wall which spatially varies in the IE- 
direction. We study one-dimensional tunneling at T = 0, 
as shown in Fig. [TJ where A and B correspond to differ- 
ent hyperfine states S z {= ±1) of the spin-1 Bose-Einstein 
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condensate. 

The main results of this paper are summarized as 
follows. We find that the domain wall structure of 
the spin-1 Bosc-Einstein condensate differs considerably 
from that in the ferromagnetic Hcisenberg model. Its 
transverse spin wave also exhibits strikingly different tun- 
neling properties. The transverse spin wave in ferromag- 
nets transmits perfectly through a domain wall [20| . In 
contrast, the transverse spin wave of the present Bose 
system is perfectly reflected by the ferromagnetic do- 
main wall in the low-energy limit. For this perfect re- 
flection, the wave function does not vanish in the low- 
energy limit. This contrasts with a conventional single 
particle in quantum mechanics, which undergoes perfect 
reflection in the low-energy limit usually due to the ab- 
sence of the wave function. The quadrupolar spin wave 
in a spin-1 Bose-Einstein condensate exhibits conven- 
tional reflection properties in the long-wavelength limit 
for which perfect reflection occurs due to the absence of 
the wave function in the long- wavelength limit. The co- 
existence of both perfect reflection and a non-vanishing 
wave function in the low-energy limit is due to the exis- 
tence of a damping mode whose damping length becomes 
infinite in the low-energy limit. 

Spin-1 Bose-Einstein condensates have another 
Nambu-Goldstone mode: the Bogoliubov excitation. In 
the low-energy limit, the Bogoliubov mode does not un- 
dergo perfect reflection when it is incident on the do- 
main wall. This is similar to the situation for a scalar 
Bose-Einstein condensate with a dark soliton. How- 
ever, energy-dependent non-monotonic tunneling prop- 
erties are observed at a finite energy. In particular, the 
perfect reflection occurs when the bound state appears 
at the domain. This differs from the dark soliton case in 
the scalar bosons, where the transmission coefficient is 
independent of energy 

This paper is organized as follows. In Sec. [Hi we 
compare the domain wall in a ferromagnetic spin-1 
Bose-Einstein condensate with that in the ferromagnetic 
Heisenberg model. In Sec. IIIH we examine the tunnel- 
ing properties of excitations through the ferromagnetic 
domain wall. 



II. COMPARISON WITH DOMAIN WALL OF 
HEISENBERG MODEL 



At T = in a three-dimensional system, the domain 
structure of a spin-1 Bose-Einstein condensate is de- 
scribed by the Gross-Pitaevskii equation [2l|, HH for the 
condensate wave function $ = (3>+i, $0, 3>-i) T (where 
the subscripts ±1,0 represent hyperfine states in the 
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FIG. 2: Calculated (a) condensate wave functions, (b) par- 
ticle density, and (c) spin density of F z as a function of x. "I>o 
is absent everywhere. £ is the healing length and is given by 
£ = h/ y/mc+p. 
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Here, we have introduced h(r, t) 



-fr 2 V 2 /(2m) 



c 0j o(r,t), h±(r,t) = h(r,t)± Cl F z (r,t) and F± = F x ±iF y , 
where m is the atomic mass, p(r,t) = $'(r, f)$(r, t) is 
the particle density, and F = $^ (r, i)S4>(r, t) is the spin 
density. The spin quantization axis of S = 1 spin matri- 
ces S = (S x , S y , S z ) is chosen to be parallel to the z-axis. 
The two coupling constants cq = 4irh 2 (ao + 2ci2)/(3m) 
and ci = 4irh 2 (a2 — ao) /(3m) are for spin- independent 
and spin-dependent interactions, respectively [22|, where 
as is the s-wave scattering length for the total spin 5 = 
or 2 channel. 

We consider a ferromagnetic spin-dependent interac- 
tion (ci < 0) and a planar ferromagnetic domain wall 
with only an x dependence (i.e., F z {x = ±00) = ±p, 
where p is the density at \x\ = 00). The condensate 
wave function for the stationary state is determined from 
((TJ by setting 3>(x,£) = e~ vt/ ' s $(x), where p is the 
chemical potential, under the boundary condition $(a; = 
-00) = (0,0,Vp) T and &(x = +00) = (yp, 0,0) T . Fig- 
ures [2(a)— (c) respectively show the spatial profile of the 
condensate wave function, the particle density, and the 
spin density of F z , which were obtained by solving the 
Gross-Pitaevskii equation with the chemical potential 
p = (co + c\)p. The result of Figs. H(a)-(c) is charac- 
terized by $0 = 0. 

The domain wall shown in Fig. [2(a) is quite different 
from that in a ferromagnet described by the Hcisenberg 
model. In the continuum limit, this spin model is de- 
scribed by 
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where the easy axis is taken to be parallel to the z- 
axis. <j = (<7 X , <r y , er 2 ) T is the spin vector, a is the 
lattice constant, and J(> 0) and K{> 0) are respec- 
tively the exchange coupling and easy axis anisotropy 
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constants. A domain wall solution to ([2]) is given by 
(<J x ,a y ,a z ) = (0, a/ cosh(x/A), a tanh(:r/A)), where A = 
\J J/K [T3, E3|- The magnitude of cr is homogeneous 
and the spin direction continuously changes around the 
magnetic wall. The present ferromagnetic domain wall 
for the spin-1 Bosc-Einstcin condensate differs markedly 
from this domain wall in the following two ways: (i) 
F XiV (oc |$o|) is absent because $o = everywhere, even 
near the domain wall; (ii) the magnitude of the spin den- 
sity is not homogeneous and it vanishes at the center of 
the wall (see Fig.^c)). 

A domain wall structure with F x>y ^ has been also 
discussed in the system we are considering [l9j . We find 
that energy of a domain wall in Ref. fl9l | with a non-zero 
value of $o is higher than that shown in Fig. [2] 

The different domain wall structures between the 
Heisenberg model and this spin-1 Bose-Einstein conden- 
sate originate from the following causes. For the spin- 
1 Bose-Einstcin condensate, the spin-density interacts 
with itself (i.e., F 2 (x)) and the spatial structure is gov- 
erned by the kinetic term of the condensate wave func- 
tion. On the other hand, for the spatial structure of 
the domain wall in the Heisenberg model, the nearest- 
neighbor spin exchange interaction (i.e., (Ver) 2 for the 
continuous approximation) in ^) is important. 

From the results in Fig. [2J the condensate wave 
function can be simply expressed by = 
(<& + i(x), 0, $_i(a;)) T , so that the Gross-Pitaevskii equa- 
tion can be simplified as 
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*±i = 0, (3) 



where c± = c ± c\. In this case, © can be regarded 
as the two-component Gross-Pitaevskii equation, if we 
regard c+ and c_ as the interaction parameters for the 
same and different species, respectively. In the ferro- 
magnetic case, we have c_ > c+, which corresponds to 
the condition for phase separation of the two-component 
Bose-Einstein condensate 24j. Figure Ufa) clearly shows 
the phase-separated domains of the two hypcrfinc states. 

Even if we add the quadratic Zeeman term to the 
Gross-Pitaevskii equation to introduce the easy axis as 
the z-axis by reference to ([2]), the situation remains the 
same. For this quadratic Zeeman effect, the absence of 
<I>o is energetically preferable [25l |26|| , That is, the ferro- 
magnetic spin-1 Bosc-Einstcin condensate has a domain 
wall solution even when no easy axis is introduced. In 
contrast, a ferromagnet in the Heisenberg model has a 
domain wall with finite thickness only when it has an 
easy axis (i.e., K ^ 0) 



III. TUNNELING OF EXCITATIONS 
THROUGH DOMAIN WALL 

The tunneling properties of excitations can be deter- 
mined by considering (fr(x,t) — exp (— i[it/h)(<5>(x) + 
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FIG. 3: Calculated transmission probability r through the 
domain wall as a function of the incident momentum when 
the transverse spin mode is incident on the domain wall. 



cf>(x,t)) (where <f> = {4>+\, 4>o> <P~i) T describes fluctua- 
tions of the condensate wave function about the mean- 
field value) and retaining terms to 0{ij>{x)) in the Gross- 
Pitaevskii equation (fT]). The ferromagnetic spin-1 Bose- 
Einstcin condensate with S z = 1 is known to have three 
kinds of collective modes [2l|, [22| . Among them, the Bo- 
goliubov mode is associated with phase fluctuations of 
the order parameter, where spin degrees of freedom is 
not crucial. The remaining two spin-wave excitations 
are associated with spin fluctuations in the S z = and 
S z = — 1 channels. While the former is called the trans- 
verse spin wave, characterized by 6F± = 8F X ± i8F y = 
v2(3>3ii0o + ^Ti^o)' the latter is called the quadrupo- 



lar spin mode, characterized by SQ^ 



2$* 



i_i and 



8Q- = 2<& + i0*L 1 ). We briefly note that the quadrupolar 
spin density is given by Q± = &S±& = 2<f>5_ 1 $ T i. 

After applying the Bogoliubov transformation, </>±i,o = 
u±\.Qe~ lEt / h — v±i e +%Et / h , we can determine the reflec- 
tion and transmission coefficients (r and r, respectively) 
by solving the Bogoliubov equation under the boundary 
condition for the tunneling problem (see Appendix [A| . 
Here, r (r) is defined as the magnitude of the ratio of the 
flux density of the reflected (transmitted) wave to that 
of the incident wave. Here, the flux density J is given by 
J = J2j=±i o Im[ujd x Uj + v*d x Vj]/m and is proportional 
to the time-averaged energy flux [27| . 

An excitation with S z = 0, corresponding to the trans- 
verse spin wave, is decoupled from other spin states 
(i.e., <f>±i), so that the transmission and reflection coef- 
ficients are obtained by solving the Bogoliubov equation 
for S z = 0. Figure [3] shows the transmission coefficient 
t of the transverse spin wave as a function of momen- 
tum. The excitation cannot pass through the domain 
wall in the low-momentum limit and the transmission 
coefficient increases monotonically with increasing mo- 
mentum. (Perfect reflection in the low-momentum limit 
is proved in Appendix [Bl) 

This perfect reflection contrasts with the case of 
the Heisenberg model ©. In the latter case, the 
eigenfunction describing the spin wave in the pres- 
ence of the domain wall is proportional to [— ikX + 
tanh(x/A)] exp (ikx) [H, H^]. That is, the incident wave 
is perfectly transmitted through the domain wall for 
all k. This phenomenon in a magnetic nanowire has 
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recently been studied by solving the Landau-Lifshitz 
Gilbert equation [20j . 

The absence and presence of the amplitude of the or- 
der parameter at the domain wall are suggestive of the 
different tunneling properties of the transverse spin wave 
for the Heisenberg model and spin-1 Bose-Einstein con- 
densates. However, a zero value of the order parameter 
at the domain wall is not related to perfect reflection. Let 
us consider the system where lim ($+1, $-1) = {\fp, 0) 

x— > ±00 

holds, but which has a region where ($+1, $-1) = (0, y/p) 
holds around x = 0. In that system, we can find points 
where <I> + i = which leads to |F| = 0. The absence 
of the spin density represents the same situation as that 
shown in Fig. [2] In this case, however, the transverse spin 
wave exhibits perfect transmission in the low-energy limit 
(as proved in Appendix |B|) . In these two cases, the ex- 
citation wave functions in the low-energy limit arc given 
by the same form of the condensate wave functions. In 
Bose-Einstein condensates, this property usually gives 
rise to perfect transmission of an excitation through a 
potential barrier in the low-energy limit [2^, [3(j. For 
the present transverse spin wave, the damping length of 
the excitation solution becomes infinite in the low-energy 
limit. This makes possible for the damping mode of the 
excitation wave function to smoothly connect with the 
zero-energy mode corresponding to the condensate wave 
function in the low-energy limit. Consequently, depend- 
ing on the configuration of the condensate wave function, 
the damping mode can remain alone in the transmission 
region where the condensate wave function is finite. See 
Appendix [Bl for more details. 

We now determine how the Bogoliubov mode and the 
quadrupolar spin mode are scattered from this domain 
wall. Hypcrfine states S z = ±1 arc coupled as in ([3]); 
however, for x = —00, the Bose-Einstein condensate of 
S z = —1 only exists, so that an excitation of S z = —1 
corresponds to the gapless Bogoliubov mode with E = 
\e(e + 2C+/9)] 1 / 2 and an excitation of S z = +1 to the 
quadrupolar spin mode with an excitation gap that has 
the spectrum E — e + 2\a\p, where e = H 2 k 2 /(2m) [22| . 
For x = 00, the Bose-Einstein condensate of S z = +1 
only exists, so that the roles for the two hyperfine states 
in the excitations are interchanged. 

Figure [1] shows the reflection and transmission coeffi- 
cients, r and t, when the Bogoliubov excitation is scat- 
tered from the domain wall. Intriguingly, reflection does 
not occur in the low-energy limit. Excitations in Bose- 
Einstein condensates are known to be perfectly transmit- 
ted through potential barriers in the low-energy limit; 
this phenomenon is referred to as anomalous tunnel- 
ing |29| ■ The Bogoliubov mode in the low-energy limit is 
a phase excitation of a Bose-Einstein condensate (i.e., a 
Nambu-Goldstone mode], which plays an essential role 
in anomalous tunneling (3J [H[. However, the absence 
of reflection in the present case differs from anomalous 
tunneling. Figure EJc) shows the matrix element of the 
local density spectral function F(x) = |<f>*w — $u| 2 [33[ 
in the low-energy limit. In contrast to anomalous tun- 
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FIG. 4: Transmission and reflection properties of excitations 
through the domain wall when the Bogoliubov mode is inci- 
dent on the domain wall. Reflection probabilities of (a) the 
Bogoliubov mode (BM) (S z = — 1) and (b) the quadrupolar 
spin mode (QSM) (S z = +1). Transmission probabilities of 
(c) the BM (S z = +1) and (d) the QSM (S z = -1). (e) A ma- 
trix element of the density spectral function in the low-energy 
regime (E = 10~ 2 c+p). (f) The real part of the wave function 
u+i, where the Bogoliubov mode is perfectly reflected. The 
coupling constant is chosen to be ci = — 3c+. 



neling, the density modes of S z = ±1 states localized 
at the domain wall appear even in the low-energy limit. 
(See Appendix [Cl for more details, where these localized 
modes are discussed.) 

With respect to scattering between a topological exci- 
tation and the Nambu-Goldstone mode, a dark soliton in 
a scalar Bose-Einstein condensate is known to be com- 
pletely transparent to the Bogoliubov mode for all en- 
ergies [3J|. In this case, the localized density mode can 
also be found in the low-energy limit, given by F (x) cx 

I, which 



tanh 2 (a;)/ cosh 4 (a;) in the dimensionless form [3 



is similar to the present case. However, in stark contrast 
to this dark soliton case, the transmission and reflec- 
tion coefficients in the present case exhibit pronounced 
non-monotonicity. In particular, the Bogoliubov mode 
is perfectly reflected at some energy points (A and B in 
Fig. IHa)). This perfect reflection is strongly related to 
the bound state of the S z = +1 state at the domain wall 
(see Fig. 0Jf)). This appears to be similar to resonance 
scattering, where the bound state increases the scatter- 
ing cross-section [35j |. However, the present bound state 
is specific to the ferromagnetic spin-1 Bose-Einstein con- 
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densate. For x < 0, (it+i,u+i) decays when the incident 
energy is lower than the energy gap of the quadrupolar 
spin mode. For x > 0, the wave functions 
behave as the Bogoliubov mode whose damping mode is 
given by exp(— kbx) for a positive energy E(> 0) (see 
also Appendix [S] and [3(| ) ■ The damping structures on 
both sides of the domain wall are different. 

For E < 2\ci\p, the plane wave solution of the 
quadrupolar spin mode does not exist and the reflection 
and transmission of the quadrupolar spin mode are both 
forbidden (Figs. 0] (b) and (d)). When the incident Bo- 
goliubov energy exceeds the energy gap of the quadrupo- 
lar spin mode (i.e., E > 2|ci|p), the transmission coeffi- 
cient of the quadrupolar spin mode (S z = —1) suddenly 
increases (see Fig.[DJd)). 

Figure [5] shows the reflection and transmission coef- 
ficients when the quadrupolar spin mode is incident on 
the domain wall. Perfect reflection occurs in the low- 
momentum limit, but excitations with finite momentum 
pass through the domain wall. At higher energies, the 
Bogoliubov mode of S z = +1 accounts for a substantial 
fraction of the transmitted flux. For perfect reflection 
in the low- momentum limit, we confirmed that the wave 
function vanishes in this limit. This is a conventional 
reflection property in the long-wavelength limit. 

The momentum transfer can be used to drive the fer- 
romagnetic domain wall. From the viewpoint of the mo- 
mentum transfer that occurs when excitations are re- 
flected, the spin wave in ferromagnets does not drive 
the ferromagnetic domain wall. This is because this spin 
wave exhibits perfect transmission that is independent 
of the energy [20(. On the other hand, in the ferromag- 
netic spin-1 Bose-Einstcin condensate, perfect reflection 
occurs for all three excitations: the transverse spin wave, 
the quadrupolar spin mode, and the Bogoliubov mode. 
The first two excitations do not play a major role in 
driving the domain wall because they exhibit perfect re- 
flection in the low-momentum limit so that they involve 
very low momentum transfer. In contrast, the Bogoli- 
ubov mode exhibits perfect reflection outside of the low- 
momentum limit (A and B in Fig. |4}. This relatively 
high momentum transfer of this mode is useful for driv- 
ing the domain wall. Spin-transfer torque can also drive 
a ferromagnetic domain wall @, Ho[. However, the for- 
malism used in this paper cannot directly demonstrate 
excitation-driven domain wall motion due to spin trans- 
fer as well as momentum transfer, because it involves 
only a small fluctuation about the condensate wave func- 
tion. The subject of future studies is to directly demon- 
strate how excitation tunneling affects domain wall mo- 
tion in the spin-1 Bose-Einstein condensate by including 
the nonlinear effect as studied in Ref. (20j . 

The simple domain wall shown in Figs. [1] and [2] can be 
mapped onto a domain wall soliton (phase separation) in 
binary mixtures of Bose-Einstein condensates [24| (see 
also (J3j) ) . In this regard, the phase separation can be ex- 
perimentally controlled by changing the scattering length 
between same or different components through the Fes- 
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FIG. 5: Transmission and reflection properties of excita- 
tions at a domain wall when the quadrupolar spin mode is 
incident on the domain wall. Reflection probabilities of (a) 
quadrupolar spin mode (QSM) (S z = +1) and (b) Bogoli- 
ubov mode (BM) (S z = — 1). Transmission probabilities of 
(c) QSM (S z = -1) and (d) BM (S z = +1). 

hbach resonance [37l [38j . The term "quadrupolar spin 
mode" is specific to the spin-1 Bose-Einstein condensate; 
however, transmission and reflection properties equiva- 
lent to those shown in Figs. 2] and [5] are easier to ob- 
tain experimentally by employing an immiscible two- 
component Bose-Einstein condensate (e.g., 85 Rb- 87 Rb 
Bose-Einstein condensates [37| and Bose-Einstein con- 
densates of 87 Rb with two internal spin states [HI). 

To create the localized Bogoliubov mode and localized 
spin- wave excitations, Bragg scattering and Raman scat- 
tering, whose light beams are shined over a razor edge, 
would be useful, respectively. After an excitation is scat- 
tered by the domain wall, if reflection occurs, excitations 
may propagate in the direction opposite to the incident 
excitation. Stern-Gerlach measurements can be used to 
determine their spin degrees of freedom. On the other 
hand, the movement and deformation of the domain wall 
after scattering can be observed by imaging the density 
because the total particle density is low at the domain 
wall (Fig.^b)). 

In this paper, to extract the essential physics, we con- 
sidered a planar domain wall and one-dimensional scat- 
tering in a three-dimensional system. Nambu-Goldstone 
modes of a Bose-Einstein condensate are known to be 
reflected and refracted in the presence of a potential 
step [HI, H(| . For junctions in Bose-Einstein condensates 
with equal densities, reflection and refraction are ab- 
sent for both the Bogoliubov mode and a transverse spin 
wave. For junctions in the (anti-)ferromagnetic Heisen- 
berg model, spin waves also transmit perfectly [4l|. It 
will be interesting to extend the present one-dimensional 
scattering problem to a three-dimensional one and to 
compare reflection and refraction properties with those 
reported in [3^, H(| and for (anti-)ferromagnets. The 
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findings in this study will be helpful for determining the 
interface physics of superfluids and scattering processes 
through domain walls that have more complex struc- 
tures than a planar wall. For scattering between low- 
lying modes and topological excitations, it is also a fu- 
ture problem to study Nambu-Goldstone mode scattering 
from a (half-)quantized vortex and a skyrmion in spinor 
Bosc-Einstcin condensates. 



IV. SUMMARY 

We studied scattering between low-lying modes and a 
topological defect in a spin-1 Bose-Einstein condensate. 
In particular, we report the structure of a ferromagnetic 
domain wall and the transmission properties of excita- 
tions through this domain wall. We found that the do- 
main wall and the tunneling properties of the transverse 
spin wave differ from those for a planar domain wall in 
the ferromagnetic Heisenberg model. We also found that 
the transmission and reflection properties are strongly 
non-monotonic when the Bogoliubov mode is scattered 
from the domain wall. 



Acknowledgments 

Authors thank D. Takahashi for useful discussions on 
zero-modes in Bose-Einstein condensates. SW thanks J. 
Ieda for discussions and pointing out references [2(], HH . 
This work was supported by Grants-in-Aid for Scientific 
Research (Grant Nos. 20500044, 21540352, 22540412, 
23104723, 23500056) from JSPS and MEXT, Japan. 



Appendix A: How to determine the transmission 
and reflection coefficients 



Considering fluctuations in the condensate wave func- 
tion Q(x, t) — cxp (—ip,t/K){Q{x) + <j>{x, t)) and retaining 
terms to 0(<fr(x)) in the Gross-Pitaevskii equation (fTJ, 
we obtain equations for determining the tunneling prop- 
erties of excitations. This appendix employs $o = 0, 
which is consistent with the system we are studying in 
this paper. For a transverse spin wave (excitation with 
S z = state) , after applying the Bogoliubov transforma- 
tion 0o = u^e^ lEt l h — VQe + * Et / h , we obtain the following 
equation 



E 



h -2ci$+ii>_ 
2ci$i 1 *l 1 -h 



(Al) 



where h = - p + c+(|$ + i| 2 + |$_i| 2 ). Ampli- 

tude transmission and reflection coefficients, To an d Rq, 
are obtained, by solving (|A1|) with the following bound- 
ary conditions: Uo = eie+ + J?o e ie~ + Aoe2e kx for 
x = — oo , and Uo = Toeie + + Boe2e~ for x = +oo, 
where hk = {2raE) 1 / 2 , Uq = (uo,vq), ei = (1,0), and 



e2 = (0, 1). Aq and Bo are unknown coefficients of the 
damping solutions, which should be determined along 
with To and Ro- 

On the other hand, as for the excitations with S z = 
±1 states, the usual Bogoliubov transformation 



u±ie 



E 



-iEt/H 



7 ±l c 



HEt/h leads tQ 



'>±1 







p 


H 











-a+ 


b+ 






fl 


-h+ 


c* 




V. 


fl 


—c 


h- 


1)1 


u 




-b*_ 


a*_ 




w 




p + 


2c+ 


|<f±i| 2 + 


c-\ 


4> 



(A2) 



where h± = ~^d 2 x - p + 2c+\&±i\- + c_|«p t i|", a± 
c + $±i$±i, b± = c_$±i$* x , and c = c_$ + i$_i. To 
determine the transmission and reflection coefficients, we 
impose the following boundary condition: 



u+i 

U_i 



^ R +1 ( q 1 ] e- ifc « x + A +1 ( 62 1 



( ° ) 



A-i j e^ x (x = -oo), 
(A3) 



u+i 

U 1 



( o ) e+lkBX + B +i ( 1 e 



+ T_i f °J e +ik * x + B-i e- KC > x (x = oo), 

(A4) 

where = (u±i,f±i), &q and kq are respectively the 
wavevectors of the quadrupolar spin mode and its damp- 
ing solution, &b and «b are respectively wavevectors 
of the Bogoliubov mode and its damping solution [36| . 
T±i and R±i are respectively the amplitude transmis- 
sion and reflection coefficients for propagating modes, 
and A±i and B±i are coefficients for damping modes, 
which should be determined in this problem. For the 
Bogoliubov mode incident to the domain wall, we set 
(u+^u'J^) = (0, a) exp (ik-Qx). On the other hand, for 
the incident quadrupolar spin mode, we set (u™ 1; u™ x ) = 
(ei,0) exp (ikqx). 



Appendix B: Proof of perfect reflection of transverse 
spin wave 

In this appendix, we prove and discuss perfect reflec- 
tion of the transverse spin wave mode in the low-energy 
limit. We will use the dimcnsionless form, where the 
wave function, energy, and length are scaled by y/p, c+p 
and £ = h/*/mc+p. Two functions S± = uq ± vo obey 
the following equation 



h ='S T T i2cilm($ + i$_i)5± = ES± (Bl) 



where h 



(±) 



h ± 2ciRc($+i$_i). Since $±i are real 



in this study, we have 



h { ±] S T = ES±. 



(B2) 
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The flux density, which is independent of x, is given by 
J cx Im[S+d x S+ + S 1 * d x S-]. For the tunneling problem, 
the boundary condition for x = — oo is given by S± = 
e ikx _|_ R Qe -ikx _j_ A Qe kx^ anc j fagfc f or ^ = -)-qq i s given by 

S± = Toe lkx ± Boe~ kx . The transmission and reflection 
coefficients r and r are now determined as r = |Xb| 2 and 
r = |i?o| 2 , respectively. 

To study the low-energy properties, we use the follow- 
ing expansion 

(T ,R ,A ,B ) =J2k n (Tt\R^,Ai n \Bi% (B3) 

n=0 

For A: | a; | <C 1, the boundary condition for x <C — 1 is 
given by S± = 1 + R { 0) ± A { 0) + k(R { Q 1] ± A^ 1 ') + - 

(0) 4. A (0)\ _l . . . ) and that for x > 1 ig given by s± = 



i^±4 0) ) 



± 4 0) + MT (1) ± B™) + M*T (0) T 4 0) ) + • • • • 
We also determine S± from the equation (|B2I) . We ex- 
pand the functions 5± as S± = ^ n=0 k n S± '. Then, 
/i^S 1 ^ ,1 = follows. Since one of the zero mode 
solutions for (|A1[) is given by (i*o,«o) = {&+i,<&-i), 



(0,1) 



= $+i ± 3>_i = </>± follow, where we took 
$±1 to be real. The other solutions are obtained as 

o(0,l) 



n \x) = 4>±{x) Jq dx'4>± 2 (x'). The behaviors for 

\x\ 3> 1 are given by (Sij , S^'P) = (l,sgn(x)) and 

(S+'n^ S -'u) = (sgn(x)7+ + x ,l- + \ x \), where j± = 
/ °° dx{<j>^ — 1). Here, we used (f>±(— x) = ±<j>±{x). 
We examine solutions linear in k, which are given by 

J±{x) — ± jD ± j+0 ± jjD ± ii + K^ ± i D ± j-|-O_|_j I ± jjj. 

By comparing the coefficients of k and x, we obtain 



(i?r^o (o) ,4 o) ,^ o) )=H, 0,0,1-0, 



(o) 
±,i 



TC± n =1 - i, 

4% =o. 



(B4) 
(B5) 
(B6) 



As a result, the excitation shows perfect reflection in the 
(0)|2 = 1). We also note that the 



low-energy limit (i.e., \R 
wave functions show lim S± = ±(1 

E->0 



(>±, that is 



lim(u ,u )/(l-i) = 



.*+!)■ 



(B7) 



has a damping solution whose damping length is 1/k. 
In the low-energy limit, this damping length becomes 
infinite and the damping mode can smoothly connect 
with the condensate wave function, which is uniform for 
|a?| 3> 1. Consequently, the propagating wave mode can 
vanish in this limit. 

The coexistence of perfect reflection and a non- 
vanishing wave function is also observed in the scalar 
Bose-Einstein condensate. In particular, we find perfect 
reflection in the low-energy limit when the Bogoliubov 
excitation is incident on the region where the condensate 
density is absent due to the presence of a potential step 
that exceeds the chemical potential [12, EH . If we regard 
the excitation in the low-energy limit as the phase mode 
of the condensate, perfect reflection in the low-energy 
limit is natural because a region with no condensate den- 
sity cannot support the phase mode. On the other hand, 
in the present case, the spin wave mode in the low-energy 
limit, which is related to the rotation symmetry of the 
spin space, can propagate on both sides of the domain 
wall, so that perfect reflection of this mode originates 
from different causes. A zero-value of the order parame- 
ter at the domain wall is specific to the present system, 
and is suggestive of perfect reflection in the low-energy 
limit. 

However, perfect reflection of this spin mode does 
not originate from the zero-value of the spin density 
at the domain wall. To demonstrate this, we assume 
that lim ($ + i,$_i) = (1,0); however, the system 

x— >±oo 

has a region where ($+1, <E>_i) = (0,1) around x = 0. 
For the tunneling problem, the boundary conditions for 
excitations are unaltered. Solutions S±'P and S^'^ 
are formally given in the same form as those derived 
above. However, they have different behaviors, which 
are given by S±'P — 1 and S±fi = sgn(x)7± + x, where 

7-t = /g°° ^a;(0± 2 — 1)- Following the same procedures 
shown above, we obtain 



/d(0) t (0) 
^0 > ± i 



(C 





(0) r (o) 
±,I' °±,II' 



|(0) d(0) 



B , 

r (l) \ 



= (0,1,0,0), 

=(i,o,*). 



(B8) 
(B9) 



The spin density fluctuation 5F, defined here as the spin 
density F = $^S$ within the first order of <fio, is given 
by SF X!y cx F z , where F z =^ 1 -$ 2 _ 1 . 

As for the usual perfect reflection in quantum mechan- 
ics, the phase of the amplitude reflection coefficient Rq 
in the low-energy limit leads to it and the amplitude of 
the wave function vanishes everywhere. However, in this 
case, the phase of the amplitude reflection coefficient Rq 
in the low-energy limit is 37r/2 and the excitation wave 
function remains in the low-energy limit when perfect re- 
flection occurs. Perfect reflection and a non-vanishing 
wave function coexist in the low-energy limit for the fol- 
lowing reason. The transverse spin mode in this system 



As a result, we have perfect transmission in the low- 
energy limit (i.e., |T , d°' l | 2 = 1), and lini£^o( u o, vq) = 
(<f>+i, $_i). The spin density fluctuation is given by 
(SF x ,SF y ) cx (F z ,0). In the present case, we can find 
points where $+i = which leads to F z = 0. How- 

ever, the absence of the amplitude of the order parameter 
at the domain walls is irrelevant to the tunneling prop- 
erties in the low-energy limit. In a scalar Bose-Einstein 
condensate, the Bogoliubov excitation can pass through 



the dark soliton [34| . This is another example in which 



perfect transmission occurs in a system where the ampli- 
tude of the order parameter is absent. 
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Appendix C: Details of Zero-Mode Tunneling 
through a Soliton 



In this appendix, we discuss effect of zero-energy 
modes on the tunneling properties by focusing on the 
phase and density fluctuations. 

First, by comparing (|A2|) for E = with ^ and its 
derivative with respect to x, we find two simple zero- 
energy mode solutions for (|A2|) . which are given by 



(«+!,«+!,«_!, »_ X ) =(*+!, (CI) 

{u +1 ,v + x,u_ u v_x) =d x {$ + i,-$* +1 ,$_ 1 ,-$*_ 1 ), 

(C2) 

where we omit the normalization factors. If we add an 
external potential term to © and (|A2[) by replacing p 
with p — V(x), the first solution f|Cl[) holds, whereas the 
second solution (|C2[) docs not. (|C2[) yields a zero mode 
solution to (|A2[) only when the system is uniform. This 
is easy to understand when we consider (spontaneously) 
broken symmetries. As for the broken gauge symmetry, 
when the condensate wave function is given by 3>(a;) = 
A(x) exp[i8(x)], the function A(x) exp {i[8(x) + d9]} ~ 
<&(x) + i&(x)d8 is another solution of Q. The factor 
proportional to dQ leads to (|C1|) . (Here, we omit the 
subscript ±1; the discussion in the present paragraph 
also holds for a scalar Bose-Einstein condensate.) Thus, 
the first solution (|C1[) is the zero mode related to the 
broken gauge symmetry and it does not contribute the 
density fluctuation. (When the condensate wave function 
$ fluctuates as $ + <j>, the density fluctuation is given by 
|$ + (j)\ 2 - |$| 2 = 2Re[$*0] + 0{(j) 2 ). In the low-energy 
limit, <f> = u — v* and the solution (jClj) is not related to 
the density mode.) On the other hand, in the absence 
of an external potential, $(x + dx) ~ §>(x) + d x &(x)dx 
is another solution of ([3]). The term proportional to dx 
leads to (|C2j) Thus, the second solution ([C2|l is the 
zero mode related to the broken translational symmetry 
due to the presence of the soliton. This gives a finite 
density fluctuation and the matrix element of the density 
spectral function is given by F(x) oc (d x p(x)) 2 where 
p(x) = |<&(x)| 2 . The profile in Fig. HJe) pertains to the 
present discussion. To summarize, when the condensate 
exhibits solitary behavior, the density fluctuation does 
appear in the soliton. 

We here discuss the absence of reflection of the Bogoli- 
ubov mode in the low-energy limit from the domain wall 
shown in Fig. @] Starting with a scalar Bose-Einstein 
condensate is instructive for understanding the charac- 
teristic properties of excitation tunneling through the do- 
main wall. In the usual manner and using the usual nota- 
tion, the condensate wave function is obtained by solving 
the Gross-Pitaevskii equation 



2m 



d 2 x + V(x)-p + g\<t>(x)f 



$(x) = 0. (C3) 




(b) 



«■) 




^ground state (V=/ 0); 
"ir-state(V=0) 
"ir-sFate(y^O) ^ 
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FIG. 6: Transmission and reflection properties of the Bo- 
goliubov excitation in a scalar Bose-Einstein condensate. We 
consider three cases: (i) the ground state in the presence of 
the barrier (V =^= 0) (solid lines in each panel), (ii) 7r-state in 
the absence of the barrier (V = 0) (dotted lines in each panel), 
and (iii) 7r-state in the presence of the barrier (V ^= 0) (dashed 
lines in each panel), (a) the condensate wave function as a 
function of x, (b) momentum-dependent transmission prob- 
ability r = |T| 2 , and (c) the matrix element of the density 
spectral function in the low-energy regime = 10 -2 . £ s is 
the healing length, which is given by £ s = h/^/rngp. For cases 
(i) and (iii), we used V(x) = 2gpexp(— x 2 /£ 2 ). 



The Bogoliubov excitation is described by 



E 



h 

gm 2 



-g<$> 2 
-h 



(C4) 



where h = -h 2 /{2m)d 2 + V(x) - p + 2g\<P{x)\ 2 . The 
chemical potential p and the boundary condition u = 
(u, v) for the tunneling problem are respectively given by 
p = gp = g\G>(x -> 0)| 2 and 

u =ae +lkBX + Rae- lkBX + A/3e KBX {x = -oo), (C5) 



u =Tae 



-\-ik&x 



B/3e- KBX {x = oo), 



(C6) 



where fee, kb, ol and f3 are given in |36| (where c + is re- 
placed by g). T and R are respectively amplitude trans- 
mission and reflection coefficients for propagating modes 
and A and B are the respective coefficients for damping 
modes, which should be determined in this problem. We 
consider three cases below. The first case is the ground 
state in the presence of the barrier V ^ 0, where the 
boundary condition of the condensate wave function is 
given by $(x — > ±oo) = + v / p. This case was studied 
as an anomalous tunneling problem (29j . The results are 
shown by the solid lines in each panel of Fig. [5] In this 
case, the Bogoliubov mode perfectly tunnels through the 
barrier in the low-energy limit (Fig. [5Jb)) and the den- 
sity mode is suppressed (Fig. El[c)). The second case is 
the 7r-state in the absence of a barrier (V = 0) where 
the boundary condition of the condensate wave function 
is given by $(x — > ±oo) = ±\ff>- This case was consid- 
ered in j3J]. The results are shown by the dotted lines 
in each panel of Fig. [5] In this case, the Bogoliubov 
mode perfectly tunnels through the soliton, independent 
of the energy (Fig.[5][b)). In contrast to the previous case, 
density fluctuation appears (Fig. [H{c)). We here exam- 
ine whether the Bogoliubov mode transmits through the 
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barrier in the low-energy limit when the system is in the 
7r-state. The results are shown by the dashed lines in 
each panel of Fig. |6] In the low-energy limit, the Bogoli- 
ubov excitation in the 7r-state perfectly tunnels through 
the barrier at the center of the soliton. In this case, the 
potential barrier suppresses the zero mode of the density. 
From these three cases, we find that the density fluctua- 
tion is irrelevant to perfect tunneling through a soliton in 
a scalar Bose-Einstein condensate. Since the presence of 
the zero mode of the phase is related to the phase coher- 
ence of the condensate, this coherence is important for 
perfect tunneling of the low-energy Bogoliubov mode in 
a scalar Bose-Einstein condensate. 

We now examine the case for the domain wall of a spin- 
1 Bose-Einstein condensate. We add an external poten- 
tial to the center of the domain wall, which suppresses 
the zero mode of the density, to investigate whether the 
perfect conversion of energy flux from S z = — 1 Bogoli- 
ubov mode to S z = +1 Bogoliubov mode occurs only 
through the zero mode of the phase. We add an exter- 
nal potential term to and (|A2[) by replacing pi with 
[i — V(x). The results are shown in Fig. [7J We find that 
perfect reflection occurs in the low-energy limit when the 



Bogoliubov mode is injected into the domain wall de- 
formed by the barrier. This result contrasts with that 
in the 7r-state of a scalar Bose-Einstein condensate. We 
confirmed that in the presence of the barrier, the matrix 
element of the density spectral function near the domain 
vanishes as the energy decreases. This result suggests 
that the localized density mode plays an important role 
in converting the energy flux of a low-energy Bogoliubov 
mode propagating in one Bose-Einstein condensate into 
that in the other condensate. 

We also confirmed that reflection is absent in the low- 
energy limit when the barrier is located in a region where 
the condensate is uniform (i.e., far from the interface be- 
tween condensates). In this case, the potential barrier 
breaks the translational invariancc so that the zero mode 
solution is not given by d x <S>±\. However, we find the 
same behavior of the matrix element of the density spec- 
tral function as shown in Fig. SJe). The breaking of the 
translational invariance by a potential barrier far from 
a domain wall is irrelevant to the appearance of the lo- 
calized density mode near the domain wall and to the 
absence of reflection. 
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FIG. 7: Transmission and reflection properties of excita- 
tions through the domain wall deformed by the potential bar- 
rier V(x) when the Bogoliubov mode is incident. Reflection 
probabilities of (a) Bogoliubov mode (BM) (S z = — 1) and 
(b) quadrupolar spin mode (QSM) (S z = +1). Transmission 
probabilities of (c) BM (S z = +1) and (d) QSM (S z = -1). 
The coupling constant is chosen as ci = — 3c+ and we used 
V(x) — 2c+pexp(— x 2 We here took the center of the 
domain wall (i.e., x — 0) as the center of the potential bar- 
rier. 



